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NOMENCLATURE 
blade section lift curve slope. per radian 
aerodynamic center 
lift deficiency function, including blade wake effects 
center of gravity 
blade section chord 
blade bending modulus of elasticity 
elastic axis 
displacement of blade k bending mode 
blade section second moment of area, flatwise bending, 
Ir 
or — 
Ip 
blade moment of inertia about flapping axis 
blade section moment of inertia about its elastic axis 
nondimensional gereralized mass of k bending mode 
blade section moment of inertia about its center of 
gravity 
nondimensional blade product of inertia about flapping 


and elastic axes 


blade bending moment at r 
blade section aerodynamic moment about elastic axis 


nondimensional aerodynamic feathering moment derivatives 


nondimensional aerodynamic derivatives associated with 


blade 


2711 





blade mass per unit length 
blade radius 
radial station along blade 
blade thrust 
velocity perpendicular to blade span in vertical 
plane at blade element 
velocity perpendicular to hlade span in horizontal 
plane at blade element 
T/R 
distance between A.C, and E.A., positive forward 
distance between C,G. and E.A., positive forward 
vertical displacement of E.A. of blade element from 
horizontal 
displacement of blade, rigid flapping mode 
Locke number pacRT 

Ip 
mode shape of blade k bending mode 
blade pitch 
nondimensional root of bending-torsion characteristic 
equation 
nondimensional rotating, natural frequency of k bending 
mode 
air density 
rotating natural frequency of k bending mode 


nonrotating natural frequency of blade torsional motion 


iv 
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SUMMARY AND INTRODUCTON 

This report presents the results of a theoretical approach to the 

analysis of the flutter of extremely flexible, helicopter rotor blades. 

The research was initiated by the fact that during experiments conducted 

at the United States Naval Academy, 1969-1970, it was found that there 

were no theoretical or experimental valves available for comparison. Some 
research had been done in the derivation of equations of motion for a 

flexible rotor blade, but actual numerical calculations had not been made. 
With this in mind, the equations of motion, incorporating three degrees of 
freedom, were derived for a rotating, flexible rotor blade. The three degrees 
of freedom used were pitching, flapping and the first elastic mode of flap- 
bending. The equations were then expanded into the characteristic equation. 
The objective of the report is to analyze the effect of the modulus of elasticity 
upon the flutter boundaries of a rotor blade. Also, at each modulus of elas- 
ticity, it is desired to study the effects of the movement of the aerodynamic 
center and the center of gravity upon the flutter boundary. 

The area of work is done entirely in the field of classical flutter. 

It is assumed that the inflow of the blade is not stalled. Hence, there is 
no wake or stall flutter to consider. 

The theoretical results consist primarily of stability plots made for 
each different modulus of elasticity. These plots are arrived at through the 
use of Routh's criteria. All of the results obtained are in very close agree- 
ment with personal observations made by the author Of experiments conducted at 
the United States Naval Academy (Reference (23)). They also agree with obser- 
Ion made by the Martin Company (Reference (12)). However, the results do 


vary somewhat from what previous research has shown to be true for fairly rigid 





blades, thus implying that more work should be done in the area of extremely 


flexible rotor blades. 
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DERIVATION OF EQUATIONS OF MOTION 


Since, in this report, it was desired to study the effects of 
modulus of elasticity on a rotating helicopter blade, it was necessary 
to derive equations of motion which would incorporate this as a 
variable. This necessitated having the equations being of more than 
two degrees of freedom in order to incorporate the bending effects 
of the blade. Three degrees of freedom were considered; the first 
torsional mode, the zeroth flapping mode, and the first elastic mode 
of flap-bending. The two flapping, modes are illustrated in Figure 23. 
Ideally, of course, an infinite number of degrees of freedom should 
be used for the analysis, but the mathematics become impossible to 
handle as can be seen in Reference (1), which only used ten degrees 
of freedom. Also, it was decided to use only three degrees of freedom 
because it was felt that these were the ones which have the primary 
effect on blade motion. 

The method of the derivation of the equations of motion vary 
from reference to reference. For example, References (2), (21), and 
(11) use LaGranges method for the derivation, whereas References (15) 
and (5) used the virtual work method. However, the method used in 
this report follows that of References (18) and (13) and is the sun- 
mation of forces or moments. 

Prior to the derivation a number of assumptions concerning the 
rotor blade were made. They are: 

(1) The rotor blade is untwisted・ 


(2) Only the first torsional mode is considered. 


(3) The rotor blade is untapered. 





(4) The rotor blade is flexible, or in other words, the first 
elastic mode is considered. 

(S) The rotor blade is of a uniform mass distribution. 

(6) The rotor blade is hinged at the root, or is articulated, 
with no hinge offset. 

Further assumptions, as needed, will be made as the derivation 
progresses. 

First, consider Figure 1, which shows the bending out of the 
plane of rotation of a flexible rotor blade. The positive direction 
for "z" will be taken as a bending of the blade upward. It is ev- 


ident that the bending moment at r due to the forces at s will be 


ward: LI Tan 2 Z reizt 5- rlds 


R > (1) 
1 تار‎ 
To arrive at the differential equation for beam bending, equation 


(1) must be differentiated twice with respect to r, which yields 


dM) - dT - mt) Zr) + 20 کت‎ 


ip a, dr” I 
(2) 


- (7) 
dr 


and from beam theory 
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Therefore, the differential equation for a blade bending out of 
the plane of rotation in a centrifugal force field is, after rearrange- 


ment 


d 


S が 
ET diz] dz) ms o*ds 
dr” due 


BC 

2 e (3) 
1 ' m aÉ + MZ = dl 
dr と 
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The form that equation (3) is in does not lend itself to a flutter 
analysis. Hence, with this in mind, a series solution, written in 


terms of normal bending modes was assumed: 
๐ ๐ 
- ٢ 
と E nu C) 9, C (4) 


For the moment, assume that the aerodynamic damping is zero -- i.e. 
dT(r)/dr = 0. Then, substituting equation (4) into equation (3) the 


bending equation becomes 


7 2 2 R 2 
> dE Pe wies 
za de 1 des y 


2 x 
+ rm dn q, * mue e zo 
dr 


assume that the motion is simple harmonic at the k rotating, undamped 


(5) 


natual bending frequency. Therefore: 
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hence ۱ 3. - X D MAN 9k 
Now, putting the aerodynamic damping back into the equation, and 
with the above assumption, the bending equation of motion becomes 
ญู 
p^ 


2 2 n " 
27 ao Vk Ik t A 0 = dl (6) 


It is now necessary to consider the effect of torsional motion 
upon the bending equation, There is an inertial term and a centrifugal 
force term due to the torsional motion which causes the coupled blade 


bending - torsion equation of motion to be 


๐ ๐ 


eo 
^. の 2 at 
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(7) 
dr 
multiplying by ny and integrating from 0 to R we get 
2 2 2 (6 2 IE 
M, 3, * Mí, (09, - (8 t 10) ), Ma, k dr (8) 
R - | 
=f Kar Ir 
where E R 2 
M, = f m dk lG ๓ ๐ 





The value of k is dependent upon the number of modes of bending 
chosen. In this report k = 2 since the rigid flapping and the first 
elastic mode are the two modes of bending being considered. 


It will be assumed that the bending motion can be described by 


9 
Z=rg+nhg (9) 
where r = rigid flapping mode 
n = first elastic mode 
By substituting (9) into (8) and nondimensionalizing the resulting 
two equations by dividing by 1,22 we get 
+ um T = 09 = Fr dT dr (10) 
As و ا‎ 
TE PN eu 
2 5 229 — อ - 12%29 Lt di dp 
E o SCAM d (11) 
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and also，the nondimensional rotating natural frequency of the rigid 
flapping mode, V,, is equal to 1.0. 

Equations (10) and (11) are the two bending equations of motion 
for a rotating, flexible rotor blade. 

Next, the blade torsional equation of motion must be derived. 
Figure 2 shows the torsional geometry used. It can be seen, from the 


figure, that the positive direction for "9" has been chosen as a pitch 





down of the blade. 


Taking moments about the elastic axis 
^ 7 x: เด R 2 
- C m(z-x6)X4m c (Tr, 8de * ( T, aede 
2 R 2 d* Í 2 dz 
+ L w, ج‎ 一 f. y mt TE kr dr - [mxo rio -x.0)dr (12) 
R 
B 
の 
substituting for z and nondimensionalizing by LA we get 


Ig + i+ le = کے‎ 0 i. 





R (13) 
= ES d Mae dr 
x2] ma 5 dk 
L R 2 R de 
where - ก ' | | ma dr + y Tod r | ED 


Equation (13), then, is the torsional equation of motion for a 
rotating, flexible rotor blade. 

It is now necessary that the aerodynamic terms of equations (10), 
(11), and (13) be expanded. From Reference (3), page 272, for a unit 


span of a thin airfoil oscillating in an incompressible flow we have 


pec [E tU 6- Ct - 2 5c)8]‏ 1- 01ھ080 
dr อ‏ 

x UE C Lu, HUTE t 5c-x) 6] = 
2 
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where 


,/"۶۶ + +0 sin Y 
: a 

Up? Z t& à R q, cos Y 

The term, C'(K), is a modified Theodorsen function (see Reference 
(16)). It has been modified to take into account the wake effect of 
the rotating rotor blade. In this report, it is assumed that C'(k) 
has a value of one. Since this report is concerned only with the 
hovering case, this assumption is a very good approximation. 

Normally, in this type of analysis, the mass terms which are pro- 
portional to 2 and D are neglected, as they shall be in this report. 
If it is desired, the effect of these terms can be included by the ap- 
propriate modification of the values of the blade inertial constants, 
but the change is usually of a negligible magnitude. 

From References (18) and (21), the aerodynamic terms can be 


expressed as 


R 
ナー (ルレ d7 gf =-m. 2 -m ۔ہہر _ غ‎ 
LOS dr : n ๓ ๑ 
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where 
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m=? CO lG G xax tasin Y (E CR) AX 
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and finally, 
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(18) 
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where 
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Although these aerodynamic expressions look very complicated, 
by assuming that the rotor is operating in the hovering state, y 
equals zero. This fact greatly simplifies the above expressions. 


02 


By assuming that the advance ratio is zero, the following 


Also, if C<<R, then MR and M^, are normally set equal to zero. 


derivatives become zero: 


Mg = 0 

Mg2 = 0 

Mg» = 0 

M = Û 
8 

Mg = 0 

M 20 
B 


Hence, the final form of the equations of motion become: the flapping 


equation 


9 2 


the flap-bending equation 


a os 
“+ m - ػ5"‎ ม เม 
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and the torsion equation 


Ck SM pp + IS, + Me 
KE CN Ly $ p + Meo 


+[M, + I +1)" Je - Ind FM, y - Le q =0 
as N 


(21) 


This concludes the derivation of the equations of motion for a 
Appendix IV gives an explanation of 


flexible, rotating rotor blade. 
the various terms in the above equations. 
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ANALYSIS OF EQUATIONS OF MOTION 


There are a couple of ways in which the equations of motion can be 
analyzed in order to obtain the flutter boundaries. By one method, the 
equations of motion have to be solved for the critical flutter frequency 
and the RPM at which the flutter instability begins. This is a very 
good approach to the problem as long as there are no more than two degrees 
of freedom. As soon as the system is increased by even only one degree 
of freedom, the mathematics involved in the solution become unmanageable, 
Without the aid of a lot of computer time. Due to the time restriction 
placed on this report this method was ruled out as just being too time 
consuming. If there is interest in this approach, References (8) and 
(19) would be of great advantage. 

A second, much simpler method of obtaining the flutter boundaries is 
to use Routh's Criteria. This, as in the first method, also requires that 
the characteristic equation be obtained. However, a solution of the equa- 
tion is not needed. According to Routh, for a system to be dynamically 
stable it is necessary that all the coefficients of the characteristic 
equation be of the same sign. Also, at the same time, it is necessary 
that Routh's discriminant also be positive. This discriminant is obtained 
from an array. The development of this array and the discriminants is 
shown in Appendix II. 

Now that the method of analysis is known, the equations of motion 
must be put in the correct form. In other words, the characteristic equa- 
tion of motion has to be obtained. Before this is done, though, it should 
be pointed out that there are no time-dependent terms in the equations of 


motion. This is because the hovering condition has been assumed and hence 


14 





the advance ratio is zero which eliminates the time-dependent terms. This 
means that all of coefficients of the characteristic equation will be 
constant, and that a conventional technique for an explicit solution of 
the equation could be used if so desired. Even though solution was not 
attempted, having constant coefficients simplified the analysis. 

To further simplify the analysis, it was assumed that the motion of 
the blade was a simple harmonic motion. Therefore, it was known that the 
harmonic perturbations of the variables 8, 0, and g would be of the 


following form: 


Vat 
B = Qe e (22) 
7م‎ 
1 (23) 
VOT 
a = Ore (24) 


Substituting (22), (23), and (24) into (19), (20), and (21), and after 


Qt 


dividing all the equations by e" , the equations of motion become: 


the flapping equation 
z z 
QN i RCo * Ro SE tlm, -Ive + Ms ๆ ๑ 7= ๐ (25) 


the flap - bending equation 
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and the torsion equation 


AS 三 Ty bo t TeV" t M8, Y 


(27) 
we | 2 z 
+ IM + + 1° [a Tr V * M,47-In$ = 0 
It can be seen that now all three equations are in terms of v, 


which is a non-dimensional root. 


y = 


IE 


If the equations were solved for V, the answer would be the critical, 
reduced frequency for flutter. Since the equation was not to be solved an- 
other parameter had to be chosen for the stability plots. This was the 
nonrotating, reduced torsional frequency, SE 

The equations can now be put into a determinant form, and the de- 
terminant expanded to arrive at the characteristic equation. This also 
shows why more than a few degrees of freedom are not very often analyzed. 
With three degrees of freedom the determinant is a 3 X 3, and the expansion 
is fairly complicated. As the degrees of freedom are increased the square 
determinant increases by one and becomes more and more difficult to expand. 
The full expansion of the determinant used in this report may be found in 
Appendix I. A quick look at Appendix I will show how unwieldy the coef- 
ficients of the characteristic equation become, with some coefficients 
having as many as thirty (30) terms. Increasing the system by only one de- 
gree of freedom nearly doubles the number of terms of the coefficients, and 
the complexity of the problem. 

After expansion of the determinant, the characteristic equation is 


found to be a sixth-order equation of the following form: 
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Once the characteristic equation was derived, it was possible to 
calculate the coefficients for the varying parameters and to establish the 
flutter boundaries. The main parameter of interest was the modulus of 
elasticity of the rotor blade. Since the equations of motion were derived 
in terms of normalized mode shapes, the only place the modulus of elas- 
ticity enters the equations is through the mode shapes. Hence, it was 
necessary to use several different mode shapes in order to study the ef- 
fect of a modulus of elasticity tending towards zero. The specific mode 
shapes used and why, is explained elsewhere in this report. 

Besides studying the effect of the modulus of elasticity, it was also 
desired to study the effects of other parameters at each respective mode 
shape. At each mode shape, two other parameters were varied. These were 
the center of gravity location and the aerodynamic center location. The 
location of both of these axes was with respect to the elastic center axis. 
The center of gravity and the aerodynamic center were varied first inde- 
pendently and then they were varied together. All of the results were then 
plotted, In all the plots the ordinate is =o which is the reduced, non- 
rotating natural frequency of blade torsional motion. By using this as a 
parameter it is possible to tell, for a particular torsional frequency, at 
what RPM the blade system becomes unstable. Also, for a particular RPM it 
1s possible to calculate the corresponding critical nonrotating torsional 
frequency. The abscissa of the stability plots was either the center of 


gravity location or the aerodynamic center location. 
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Due to the complexity of the coefficients of the characteristic 
equation, all the numerical calculations were done on a computer. 
The computer used was IBM 360 Model 91. The computer was located at 
the Princeton University Computer Center at Princeton University. Since 
the program used was just a simple mathematical routine, it is not in- 
cluded in this report. All of the results were transferred to graphs, 


which are included, however. 
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MODE SHAPES 

As has already been seen, the equations of motion of the rotating 
rotor blade were derived in terms of normalized mode shapes. It should 
also be noticed that the modulus of elasticity does not appear explicitly 
in the equation, but rather is incorporated into the mode shapes. This 
made it necessary to choose various mode shapes which would be representa- 
tive of the varying modulus of elasticity. It was possible to obtain ex- 
act mode shapes for two limiting cases (see References (18) and (14)). 

The two limiting cases are 

(1) A rotor blade with the modulus of elasticity equal to zero. 

(2) A rotor blade that is considered to be a typical rigid rotor blade. 

It was not possible to find exact mode shapes for a modulus of elas- 
ticity between case (1) and case (2). However, by graphing the above 
two mode shapes and then drawing curves which fell between these two cases, 
it was possible to obtain other mode shapes. See Figure 3. Only 
two intermediary mode were chosen, making a total of four (4). 

The reason that an infinitely stiff blade was not chosen as the lim- 
iting case was that it is nearly impossible to achieve such a situation 
with a long, thin rotor blade. Therefore, the mode shape of a typical 


rotor blade was used. Mathematically the mode shape can be expressed as: 


2 
= 3X- 4X (29) 
This mode shape was then substituted into the nondimensional derivatives 


giving the following: 


= Me- Lf -+#*| SEN 
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These derivatives were then used to calculate the coefficients of the 
characteristic equation. 

Now look at case (1); the rotor blade with the modulus of elas- 
ticity equal to zero. The mode shape for this case was evaluated in a 
Ph, D. thesis written at Rensselaer Polytechnic Institute. (Reference 


(14)). It was found that the mode shape was of the following form: 


3 
KT [3x- 5X7] GO 
With this as the mode shape the non-dimensional derivatives which are 


affected by the mode shape now have the following values: 


The rest of the derivatives do not change in value. 
After much drawing of curves and curve fitting, it was possible to 
arrive at two other mode shapes which fell between the two limiting cases. 


The first of these intermediary shape was of the form: 


(31) 
k= 3x- 5x? 
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with the corresponding derivatives having the values: 


-X 
حرط‎ 24 
y 
RE LL 
3% 
Moal” 24 


The second intermediary shape was of the forn: 


EE Tax A (32) 


R = 
and its corresponding derivatives had the values 
y» (E 
MO 
_49S 
74 7 4 80 
-AY 
Mer fo 


One interesting fact discovered while arriving at the mode shapes 
was that WI the nondimensional rotating natural frequency of the second 
bending mode, did not change as the modulus of elasticity changed. (see 
References (14) and (20)). Hence, it was possible to use the same fre- 
quency throughout all the calculations. In this case the frequency turned 
mut to be 2.5. 

The reason there is no change in the frequency is that it is more 
strongly dependent upon mass distribution than it is modulus of elasticity. 
And since all the cases in this report have the same mass distribution, 
i.e. uniform, the frequency of the second bending mode, or in other words, 
the first elastic mode, remains constant. 

It should also be mentioned that the normalized mode shapes are also 


strongly dependent upon mass distribution. Although this fact was not 


considered in this report, it would be worthwhile to do so. 
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RESULTS 
‘Overall Results 

As we previously mentioned, the coefficients of the characteristic 
equation and the values of Routh's discriminants were all calculated; and 
then, using Routh's criteria, the flutter boundaries were plotted. For 
all the calculations, one rotor blade type was used. The rotor blade used 
was a model rotor blade developed by Cornell Aeronautical Laboratories for 
some of their experiments. (see References (6) and (5)). The represent- 
ative data for this model blade is given in Appendix III. 

For each mode shape, four different situations were analyzed in order 
to provide a complete flutter analysis of the rotor blade for the para- 
meters involved. The cases investigated were: 

(1) The aerodynamic center and the center of gravity were coincident 
with the elastic axis. 

(2) The aerodynamic center was coincident with the elastic axis, and 
the center of gravity was varied in location forward and aft of the elastic 
axis. 

(3) The center of gravity was coincident with the elastic axis, and 
the aerodynamic center was varied in location forward and aft of,the elastic 
axis. 

(4) And finally, both the center of gravity and the aerodynamic center 
were varied in location forward and aft of the elastic axis. 

Each case shall be discussed separately, and then, an overall analysis 
shall be given. 

First, let's consider case (1), where all three axes are coincident. 


It was not possible to plot a graph of the dynamic stability of this system 
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in a manner which could easily he understood. Therefore, the results 
were put in table form and are shown in Table 1. In the table, the mode 
shapes are numbered 1 through 4, with mode shape number 1 being the most 
rigid case and mode shape number 4 being the most flexible case. The 
first thing that will be noticed is that as the modulus of elasticity is 
decreased (the blade becomes more flexible) the RPM at which the rotor 
blade becomes unstable is decreased (an increase in I is a decrease 
in RPM). 

From Reference (22) it was found thata typical operating value for 
コグ 7 is 8.0. Hence, it immediately becomes evident, that with an ex- 
tremely flexible rotor blade the blade would always be just on the verge 
of instability, requiring only a slight abnormal operating condition to 
cause an unstable flutter. To make matters even worse, there were some 
observations made by the Martin Company, of extremely flexible blades 
(Reference (12)) which showed that as the modulus of elasticity of a 
rotor blade is decreased, the nonrotating natural frequency of the blade 
torsional motion, uo, also decreases. This means that will also be 
reduced and hence, the stability margin is decreased. All of this points 
to the problems designers have with extremely flexible rotor blades. 

In the second case, the aerodynamic center and the elastic axis were 
coincident. The results are shown in Figure 4. Once again it was ob- 
served that as the modulus of elasticity was decreased, the region of in- 
stability was increased. It can be seen that there was a gradual increase 
in instability for all mode shapes as the center of gravity was displaced 
further and further from the elastic axis. However, when the center of 


gravity was located more than 20% chord from the elastic axis the system 
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became unstable for all values of =~ analyzed. Other References 


do not show this marked instability due to the fact that in the majority 
of the analyses the center of gravity is not moved as far as 20% chord 
away from the elastic axis. The reason for this sharp instability is 
discussed elsewhere in this report. 

In case (3), the center of gravity was coincident with the elastic 
axis. This says that the rotor blade was mass balanced. Again it was 
observed that a reduction in the stiffness of a rotor blade increases the 
region of instability. It was also found that as the aerodynamic center 
was moved aft of the elastic axis, the flutter boundary began to fall off, 
whereas if it was moved forward of ‘the elastic axis a plateau for the flutter 
boundary was reached. See Figure 5. This would suggest that it is pre- 
ferable to have the aerodynamic center aft of the elastic axis. That is 
the same conclusion reached by Goldman in Reference (12). Contrary to 
case (2), and as will presently be seen, case (4), case (3) does not have 
the sharp upward break in the flutter boundary. As might be suspected, 
this was because the rotor blade was mass balanced. 

Figures 6-19 show the final results of case (4), in which both the 
center of gravity and the aerodynamic center were varied. As was expected, 
it was found that in this case, also, that reducing the stiffness of a 
blade increased the instability. It was further found, again, that it was 
preferable to have the aerodynamic center located aft of the elastic axis. 
(see Figure 20). This was the same conclusion reached when case (3) was 
analyzed. Case (4) also illustrated the same strange break in the flutter 
boundary that was found in case (2). For some reason there appeared to be 


a limiting center of gravity position, beyond which the system.was unstable 
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Since an optimum positon for the aerodynamic center had been found, 


for all values of 


it was deemed desirable to also find an optimum position for the center of 
gravity. With this in mind, Figure 21 was drawn. From this figure it can 
be seen that the region of instability will be the smallest when the center 
of gravity is coincident with the elastic axis. By looking at the equa- 
tions of motion it can be seen that in this situation the torsional flutter 
is uncoupled from the flapping and flap-bending flutter. (see Reference 
(19)). As soon as the rotor blade is no longer mass balanced the pitching 
and the flapping of the rotor blade become coupled and the instability is 
increased. It was also observed that a forward center of gravity move- 
ment was preferred over a rearward movement. This fact is further supported 
by observations made by Goldman in Reference (12). 

In no other references is there any mention of forward center of 
gravity movement made. And if the rotor blade is fairly rigid there is no 
real reason to consider this fact. However, if a rotor blade being designed 
is extremely flexible, the center of gravity location could be the source of 
many problems. 

From personal observations made by the author of some experiments con- 
ducted at the United States Naval Academy with extremely flexible blades, 
Reference (23), the critical importance ofcenter of gravity location was seen. 
The rotor blades used had originally been constructed with the idea of having 
the center of gravity and the elastic axis coincident. However, it was found 
that actually the center of gravity was aft of the elastic axis. At the 
time this was not of too much concern. However, upon conducting the experi- 


ments it was not possible to keep the blade stable. This theoretical analysis 
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From personal observations made by the author of some experiments con- 
ducted at the United States Naval Academy with extremely flexible blades, 
Reference (23), the critical importance ofcenter of gravity location was seen. 
The rotor blades used had originally been constructed with the idea of having 
the center of gravity and the elastic axis coincident. However, it was found 
that actually the center of gravity was aft of the elastic axis. At the 
time this was not of too much concern. However, upon conducting the experi- 


ments it was not possible to keep the blade stable. This theoretical analysis 
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now explains why this was so. For flexible, i.e. a very small modulus 
of elasticity, the center of gravity should be coincident with or for- 
ward of the elastic axis. 

All of the above results can be reduced down to a few certain facts 
that must be kept in mind when designing an extremely flexible rotor blade. 
Number one is that it is desirable to have the aerodynamic center aft of 
the elastic axis. Secondly, it was found that the center of gravity lo- 
cation is critical and should be conicident with or forward of the elastic 
axis. It was also found that as the modulus of elasticity is reduced, the 
nonrotating natural frequency of blade torsional motion, Wo» becomes an 
important factor. The reduction of wo that occurs with a reduction of 
modulus of elasticity tends to place the operating mode of the rotor blade 
a lot closer to instability. This also explains why it is desirable to 
have the center of gravity forward of the elastic axis. The extra margin 
of stability that is gained is needed. 

One other observation made in this analysis concerned the effect of 
center of gravity displacement from the elastic axis upon the torsional 
frequency. From the figures it can easily be seen that as the center of 
gravity is displaced further and further from the elastic axis, in either 
direction, the torsional frequency required for flutter is increased. 
Flutter Phenomena 

As was found in discussing the overall results, case (2) and case (4) 
of the analysis disclosed a very sharp upward break in the flutter boundary 
that occured at a center of gravity location of more than 20% chord dis- 
placement from the elastic axis. Since this break in the flutter curve 


occured within one chord length it was considered necessary to analyze it 
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and see if the reason for its occurrence could be discovered. First, it 
was important to find out whether it was the coefficients of the charact- 
eristic equation or one of Routh's discriminants which was displaying the 
instability. 

After going back over all of the data, it was found that in all 
instances it was the coefficients of the characteristic equation which were 
displaying the instability. Also, it was observed that it was only the 
coefficients of the two highest order terms, or in other words the "A" co- 
efficient and the "B" coefficient if the equation is written in the form of 
Equation (28). The next step in the analysis was then to find out which 
term of the coefficient caused the instability, or speaking mathematically, 
which term caused the coefficient to become negative. 


The two coefficients may be expressed as follows: 


A= II, -I -I I (33) 
B= 2 Dx 15272 + Ly Ma "D + Ly Ma + 2 Mg 
+ Im, اح‎ (34) 

Although all cases were analyzed, the case that will be discussed in 
this report will be the case where the modulus of elasticity equals zero. 
The results are the same for the rest of the cases, but with a little dif- 
ferent numbers for the coefficients. 

Using the derivatives for the case of zero modulus of elasticity, and 
substituting in the corresponding numbers, the coefficients may now be 


written in the following form: 


A= ,000287 - .0156X_ 7 ‚08444, (35) 
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c 
B= -,013 Y, t ,o£3 X, XA- , 113YXp t 0001 - ,015 X, 
+ 0694," + .000386-,043XÍ + ‚000087 (36) 


Let's consider the "A" coefficient first. A look at the numerical 
equation shows that it will be due to the second and third terms that the 
equation first becomes negative and hence the system becomes unstable. 
These two terms, which are dependent uponcenter of gravity position, are 
E and I ٭‎ Both of these terms are the nondimensional blade product of 
inertia about the flapping and elastic hinges. Going back to the original 
equations of motion, it was found that it was the inertial moments about 
the pitching and flapping axes which gave rise tothe above terms. These 
terms are also called the mass unbalance terms since if the rotor blade is 
mass balanced they are zero. 

Now let's look at the "B" coefficient. It can be seen that this co- 
efficient is dependent upon both center of gravity position and aero- 
dynamic center position. However, it was found that the only time that the 
"B" coefficient was important was when the aerodynamic center was aft of 
the elastic axis and the center of gravity forward. This meant that the 
secomd term was the most important term. It was also found that I. had the 
primary effect, as it did in the "A" coefficient. 

As can be seen, in all cases the sharp break in the flutter curve was 
due to the mass unbalance terms, I, and I>: From looking at the equations 
of motion it can be seen that the instability occurs due to a coupling of 
the inertias between the pitching and flapping of the rotor blade. Also, 
it should be mentioned that these two inertia moments arise due to either 
a pitching acceleration or a flapping acceleration. 


That, then, explains the reason: for the sharp upward break in the 
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All this shows that a divergent instability becomes of no more 
importance for an extremely flexible blade than it was for a fairly 


rigid rotor blade. 
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CONCLUS LONS 

Based on the theoretical results previously discussed, the following 
basic conclusions were reached. 

1. As the modulus of elasticity of a rotor blade is decreased the 
region of instability increases. Since, at the same time, 7 is decreasing 
there can be a very real problem with thedynamic stability of an extremely 
flexible rotor blade. 

2. In the design of an extremely flexible rotor blade, it is desir- 
able to have the aerodynamic center located aft of the elastic axis in order 
to obtain dynamic stability. 

3. Also, when designing an extremely flexible rotor blade it is 
necessary to have the center of gravity coincident with the elastic axis 
or forward of it. This is contrary to classical flutter analysis, but is 
a necessary condition if the blade is to be extremely flexible, due to the 
decrease in 0 ・ 

4。 Coupling of the mass unbalance terms, or in other words, the iner- 
tial effects due to pitching and flapping acceleration causes a sharp up- 
ward break in the flutter boundary if the center of gravity is displaced 
more than 20% chord away from the elastic axis. This occurs whether the 
center of gravity is forward or aft of the elastic axis. 

5. A change in the modulus of elasticity of a rotor blade has no 


effect on the divergence boundary of the blade. 
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RECOMMENDAT TONS 

Due to the time restriction on this report there were many areas of 
the flutter problem that the author was not able to look into. With this 
in mind, the following are several recommendations for future study in 
tnis field. 

l. Since this report was strictly theoretical, an experimental 
analysis should be made to verify the results obtained herein. 

2. It would be informative to do both theoretical and experimental 
research on the effect of extremely flexible blades upon forward flight. 

3. Although this report only considered classical futter, the other 
two types of flutter - wake flutter and stall flutter - should certainly 
be looked at. 

4. One thing that probably will be a big factor in the construction 
of extremely flexible rotor blades is a tip mass. Therefore, the effect 
of a tip mass upon the flutter boundaries of a rotor blade is in itself a 
big subject to consider. Also, a mass distribution, other than uniform 


might be considered. 
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MODE SHAPE 





TABLE 1. DYNAMIC STABILITY OF CASE (1) 


DYNAMIC STABILITY 


unstable for “%/N < 5, stable elsewhere 
unstable for O < 6, stable elsewhere 
unstable for ۳۲ < 6, stable elsewhere 


eo 


unstable for Q < 7, stable elsewhere 
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FIGURE 2 BLADE TORSIONAL 
GEOMETRY 
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memo + LxI2M; -Mo Mà t nm Ma + Me må t C lm 
2 ง 1 ท t Ix more + Lira Te me 


t [2 Meme ~L Lem, e IL) D M sma mo peas má 
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+ LM; I TO "Me md 
F=-2T; DT, Ve + mon, M; me a no ne mq NA; 
= ma Mz - Ty lat I me M; ۲۰۰ ٠ت‎ 5.۱۸۰۱۷۷۰۲۶ 21 12 
RES (ETA mor my Mg * Lx Ma ma t Mo må mg f Ln vm 
+I (E) m mi m. M ni Vv dee to laine. t+ Lo Mo + جا‎ 
o TL GI + Me ma tMi morma ~ LA Ma me -Momi -Lm 
0 0 
F= - Ly2 mm, + Lx mo ma - Llamo Ms Ve tO Le La ^m, - Ir 
Us E IM; Y 1 ye "mei me 一 Luz Lam, + RE Me “ne Ve t ک0‎ E 
+ TJ, (4) m VE me. Ms + Ixa M; t Momi + Lang E má 
+l Me Vz - Jx mer ma 
G- Tlr 1 nee, HT, NE ۹۹۶۶ 
000 vr 
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APPENDIX II--ROUTH'S ARRAY 
To establish stability of the rotating helicopter blade, with 
three modes of motion, Routh's criteria was used. In order to use 
Reuth's Criteria, Routh's array had to be set up. The form of the 


characteristic equation was: 


AD مرک‎ DO t EV FF? 1 G-0 (111) 
This leads to an array of the following forn: 
ve A C E G 
y> B D ۲ 
vi A, B, G 
v? A, B, 
oF Az G 
E A, 
ง ป G 
where 


A, = BC - AD 


1 
8 

B, - BE - AF 
B 

A, = AD - BB, 
Ay 

B, = A,F - BG 
A 

An = AB) - AB, 

A, 
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and 
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APPENDIX III--ROTOR BLADE DATA 
Cornell Aero Lab Blade #1 
NACA 0012 Airfoil 
Radius - 48 inches 


Chord 3.5 inches 


Vere e 2.5 pounds 


= 7000479 slug feet e 
I, - 1.66 slug feet : 
M - 2.46 
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APPENDIX IV--TERMS OF EQUATIONS OF MOTION 


I, Flapping and Flap-Bending Equations 

1. B and g terms. 

These terms arise from the z - component of the centrifugal force 
at each blade element. A virtual deflection of QB or Og results in 


energy being expended which results in a force for the whole blade of 


a I, 6 (TV1) 
M, aay) “y (1V2) 
2. 0 term. 


Due to a twist of the blade, the centrifugal force of a blade element, 


Mrü^dr, will act with a moment arm, - «Arê, ۲۷۰٠٣٦٢71۶677 7 ٤ 


j 2 
|j mod, AX d 


3. 8 and g term. 
At each blade element there is an intertial force of 
mredr 
ہہ‎ ۷۳۴ 
which results in a force of 


ES Ze (IV4) 
l6 
M, $ (IVS) 
4. 9 term. 


Duce oO a pitching acceleration there is an inertial effect which re- 


sults in the following force. 


R T 
= MEN, Xp dr e (IV6) 
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5. 6 term, aerodynamic. 
This term is due to a twist angle, 8, which causes energy to be 


transmitted by the airstream to the blade giving a force of 


2 (LV7) 
I, e eae 8 


and 


È 
0+ e (V8) 


6. B and g terms. 


The flapping velocity at each blade element causes a relative change 


in the angle of attack which results in the following forces: 


We (TV9) 
NN 
am 4 (1V10) 
ar (IV11 ) 
Man, P (IV12) 


RI. Sprtehing Equation 
938 tern. 
Due to a virtual twist of the blade, the z - component of the cen- 


trifugal force of a blade element causes a force of 


(IV13) 


R 
2 
e Ar qr 
7 nn, Xp dr q 
2. 0 terms. 
This term arises from the "tennis racket" effect. Each partical of 
mass, dm, of a blade section is acted upon - due to a twist deflection - by 


a component of the centrifugal force in the plane of the section. The re- 


sulting force is 


(ro GE Al 9 (IV14) 
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5 ง B term. 
Due to a Virtual deflection of a hlade section, the inertia force, 


mrBdr , results in a force of 


mat xc dr € (IV15) 


Viwseterce arises from the inertia force due to a twist of the blade. 

- 2 

I,f2 e (IV16) 
Sr B term. 


1 ท ง “ๆ term xs the result of a couple produced by the 


pitching acceleration of the rotor blade. 


۷ ENE (1V17) 


~ 


CE O term: 
EE tern is tne result of the energy transmitted to the 


blade by the lift forces moving through a virtual displacement. 


Ns E, ก (LV18) 
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